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We describe an explicitly discrete model of the process of DNA branch migration. The model matches 
the existing data well, but we find that branch migration along long strands of DNA (N > 40 bp) is also 
well modeled by continuum diffusion. The discrete model is still useful for guiding future experiments. 
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Lattice methods have proved to be very use- 
ful for studying particle interactions. As reported 
at this meeting, lattice gauge theory calculations 
now seem to be able to predict the mass of the 
"glueball," and are leading to a better under- 
standing of both chiral symmetry breaking and 
quark confinement. It is implicit in such calcula- 
tions (and sometimes explicit) that one is to take 
the continuum limit, but one can, of course, ap- 
ply lattice field theory methods to systems which 
are inherently discrete, without taking the con- 
tinuum limit. DNA branch migration is a simple 
example. 

DNA (deoxyribonucleic acid) transmits ge- 
netic information when the molecule is replicated 
and transmitted to daughter cells. Rare errors in 
replication give rise to mutations. (The adapt- 
ability of organisms with these mutations gives 
rise to evolution.) Mutations do not stay in 
the DNA segment (chromosome) in which they 
arose; they can transfer to other chromosomes by 
recombination. In recombination, two DNA seg- 
ments encoding the same genes effectively break 
at homologous points, and then rejoin or recom- 
bine by fusing partner DNA segments. When 
viewed in more detail, one finds an intermediate 
state (a "Holliday structure" ) in which two DNAs 
are joined by the. e xchang e of only one strand 
from each duplexEl (|Fig. l b). The branch point 
which joins the two helices is mobile because of 
t he seq uence similarity of the homologous DNAs 
( Fig. 1 :) . A better understanding of branch mi- 



gration will lead to a better understanding of 
DNA recombination and the propagation of ge- 
netic variation which drives evolution. 



Laboratory experiments can now replicate 
DNA branch migration through short defined 
segments of DNA with controlled initial condi- 
tions.crB Our goal is to develop a model of these 
experiments which preserves the discrete nature 
of DNA. We focus on the motion of the branch 
point and ignore the helical structure of DNA. (A 
related model has been described by Fujitani and 
KobayashiJj) This can be thought of as a series 
of linked chemical reactions in a "compartment 
model. "H Imagine that all of the molecules with 
branch points in a particular position are kept 
in one compartment, and that when the branch 
point on a given molecule moves to a neigh- 
boring position that molecule is removed from 
its compartment and added to the appropriate 
neighboring compartment. If there are N pos- 
sible positions for the branch then there will be 
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(a) (b) (c) 

Figure 1: Schematic of single strand exchange in 
DNA to form the Holliday structure (b) , followed 
by branch migration (c). 
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N + 1 compartments, the last being for molecules 
which have separated. Ignoring this last com- 
partment, let the n-th component of the vector 
C(t) represent the concentration of molecules 
having the branch in the n-th position at time t, 
and let (3 represent the rate of transport in the 
forward or backward directions (assuming these 
are equal — it is easy to generalize to the case 
where they are not). The dynamics of this sys- 
tem is then described by a set of N coupled first 
order linear differential equations, which can be 
written in matrix form as 



dC_ 
~dt 



(3a 2 MC , 



(1) 



where the N x N matrix M is 



M 



/-I 
1 



V 



-2 1 
1 -2/ 



(2) 



All blank entries are zero, and we have introduced 
a as the spacing between successive positions 
along the DNA. 



The similarity of Eq. (2) to the lattice second 
derivative is obvious, and so it should come as no 
surprise that the continuum limit of this model 
is the diffusion equation. However, we seek 



solutions to Eq. (1) before the continuum limit is 
taken. The procedure for finding these solutions 
is similar to the continuum case. First, one 
can remove a time-dependent factor by writing 



C{t) 



-t/r 



v, where v is an eigenvector of M. If 



the n-th component of v is written as v{n), then 
the eigenvectors of M have the form 



v (n) = cos (kna + <f>) , 



(3) 



where n = 0, 1, . . . , N — 1, and where k and <f> are 
yet to be determined constants. The eigenvalues 
are 



-4 9 / ka 



from which it follows that 
-1 



(3a 2 A 



(4) 



(5) 



In the continuum case the undetermined con- 
stants k and <f> would be determined by boundary 
conditions, but there are no "boundaries" to a 
matrix. Nevertheless, it is useful to imagine that 
there are extra compartments at either end of the 
system. Eq. (3) is only an eigenvector because 
it fits the +1, -2, +1 differencing pattern along 
the diagonals of M. This pattern must be main- 
tained for the eigenvalue equation to be valid. To 
do so at the upper left corner of M we imagine 
an extra point at n = —1 such that v{— 1) = v(Q), 
which forces <fi — ±ka. This says that the lattice 
derivative is zero, so in analogy to the continuum 
case we call this a "lattice Neumann" boundary 
condition. To maintain the differencing pattern 
at the lower right corner of M we imagine an ex- 
tra point at n = N, such that v(N) — (this 
is the final compartment that we ignored ear- 
lier). This "lattice Dirichlet" boundary condition 
forces the unknown constant k to be 



(2rn+ I) 7T 
(27V + 1) a 



(6) 



where m = 0, 1, 2 . . . N — 1. There are exactly N 
distinct eigenvalues, which we will label v n 

is 



The most general solution of Eq. (1 



linear combination of all of the solutions, 



C{t) 



JV-l 
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m=0 



-t/r„ 



(7) 
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Figure 2: Data from an experiment measuring 
the fraction of branched DNA converted to sepa- 
rated final product, along with the best fit curve 
from the model, [Eq. (9)|, for N = 40. 



where the (Xyyi Eire constants, and where r m is re- 
lated to the m-th eigenvalue via Eq. (5) . The 
values of the Q"ui 3X6 determined from the ini- 
tial conditions, using the fact that the different 
eigenvectors v m are orthogonal. In the exper- 
iments we are modelings the initial conditions 
are such that all molecules begin with the branch 
point at one position, adjacent to a reflective 
barrier created by a divergence in the DNA se- 
quences. These initial conditions are represented 
by C(0) = Co<5n,o- Substituting this into Eq. (7), 
taking the dot product with v m i and using the 
orthogonality of the v m yields 



4Cn 



2N + 1 



■ cos 



ir (2m + 1) 
2 (2N + 1) 



(8) 



The experiments we model only measure the 
concentration of separated molecules at a given 
time, not the number of molecules having the 
branch point in a given position. To compute 
this observable we take the initial concentration, 
Co, and subtract from it the sum over all of the 
branched molecules left in the system. This gives 



F N (*) = C 



2N- 



N-l 
m=0 



(-1/ 



-t/l 



cos- 



sin 



/ 7r 2m+l\ 
^ 2 2JV+1 J 



2m+l 
2JV+1 



(9) 



Aside from the number of steps, N, there are 
only two free parameters in this expression, Co 
and (3 (wh ich is related to the T m by Eq. (5) ). 

FIgl 2 shows measurements of the concentra- 
tion of final product (separated molecules) from 
an experiment using a DNA junction in which 
the branch point must migrate 40 base pairs (bp) 
in one direction before coming to the end of the 
strands. The errors include systematic errors, 
not just statistical variations. The curve in the 
figure is a non-linear least-squares fit to [Eq. (9)| 
with N = 40. The fit is quite good. One notable 
characteristic of branch migration is the "lag" at 
the beginning of the experiment — no material is 
produced as final product for several minutes af- 
ter the experiment has begun. The theoretical 
curve reproduces this behavior well. 
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Figure 3: Difference between the best fit discrete 
and continuum measures of the final product. 



Is it possible to see the discrete nature of the 
motion of the branch point? In particular, can 
we tell the number of base pairs traversed in a 
single step of branch migratio n? Un fortunately, 
the answer is no. The data in Fig. 1 are also fit 



well by Eq. (9) with TV anywher e betwe en 4 and 
40, as well as by the analog ue of Eq. (9) obtained 
from continuum diffusion. Fig. 3 show s the dif- 
ference between Fjy(t) from Eq. (9) and F(t) 



obtained from continuum diffusion (assuming the 
same Co). Except for very short molecules, the 
differences are smaller than the size of present 
experimental errors. Nevertheless, Fig. 3| shows 
that the discrete nature of branch migration will 
be most noticeable within times less than To. 
This is useful knowledge for designing future ex- 
periments. We therefore expect that our discrete 
model will lead to a better understanding of the 
underlying dynamics of branch migration. 
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